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Abstract 

We formulate a quantum formalism for the statistical mechanical models of discretized field 
theories on lattices and then show that the discrete version of cj> 4 theory on 2D square lattice is 
complete in the sense that the partition function of any other discretized scalar field theory on an 
arbitrary lattice with arbitrary interactions can be realized as a special case of the partition function 
of this model. To achieve this, we extend the recently proposed quantum formalism for the Ising 
model |H and its completeness property (2| to the continuous variable case. 

1 Introduction: 

The understanding that a single partition function can describe different phases of matter, is rather re- 
cent and indeed as late as 1930's, there was not a consensus among physicists that a partition function 
can give a sharp phase transition. The works of Kramers, Wannier, Onsager [|3]|4]|5| and others gradu- 
ally established beyond doubt that in the thermodynamic limit, singular behavior and phase transition 
can arise from a single partition function based on a single model Hamiltonian. For example the Ising 
Hamiltonian can describe both the ordered phase of a ferromagnet and the disordered phase of a para- 
magnet. Near the point of second order phase transition, even the details of the model Hamiltonian do 
not matter and only some general properties like symmetries are important. 

Decades of works on statistical mechanics models, inspired by the above general understanding, 
has revealed even further fruitful relations between different models. Besides the well-known duality 
relations between the low and high temperature phases of the Ising model, one can also mention the 
so-called vertex models, which reduce to other models in different limits. 

One can ask if there are certain statistical mechanical models which are complete, in the sense that 
their partition function reduce to the partition function of any other model in a suitable limit? If this 
turns out to be the case, then we can imagine a very large space of coupling constants and one single 
Hamiltonian, i.e. the Ising model with inhomogeneous couplings, so that when we move through this 
space, we meet new phases and new models which at present are thought to be completely unrelated. 
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This will then be another forward in the unification program mentioned above. 

It seems that the answer to the above question may be positive. Recent results |[Tll2ll6l [81171 l9l[l0l 
[TT1 IT2l [T3ll brought about by merging of ideas from statistical mechanics and quantum information 
theory, give positive clues in favor of the above idea. These investigations (HE] have been made pos- 
sible by establishing a link between statistical mechanics and quantum information on the one-hand 
and the new paradigm of Measurement-based Quantum Computation (MQC) and the universality of 
cluster states for MQC on the other SH [Bl [161 [HI [Til [T9l |20l IHl |22l lil . 

In a series of recent works, it has been shown that the Ising model on two dimensional square 
lattices with complex inhomogeneous nearest-neighbor interactions, is complete in the sense that the 
partition function of all other discrete models with general fc-body interactions on arbitrary lattices 
can be realized as special cases of the partition function of Ising model on a square lattice which is 
polynomially or exponentially larger than the original lattice. The starting point of these developments 
was the observation in fT| that the partition function of any given discrete model can be written as a 
scalar product, 

Z G {J) = (a|* G >, (1) 

where (a\ is a product state encoding all the coupling constants J and \^f G ) is an entangled graph 
state, defined on the vertices and edges of a graph, and encoding the geometry of the lattice. The core 
concept of completeness is the fact that the 2D cluster state is universal. More concretely we know 
that the graph state \Q G ) corresponding to a graph G can be obtained from an appropriate cluster 
state l^n) corresponding to a rectangular lattice (denoted by □), through a set of adaptive single- 
qubit measurements M. The measurements being single-qubit can be formally written as |q;m)(qm|, 
where |«m) is a product state encoding the qubits the bases and the results which have been measured. 
Thus the totality of measurements Ai transforms the cluster state as follows 

!*□) — > \a M )(a M \^n) = \um) ® |*g), (2) 

which shows that after disregarding the states of the measured qubits \o.m), what is left is an appro- 
priate graph state 

|* G ) = (a u \* n ). (3) 

Combination of this relation with (fl}, leads to the completeness result mentioned above. That is one 
writes 

Z G (J) = (a\* G ) = (a,a M \* a ), (4) 

and notes that (a, au\ now encodes a set of generally inhomogeneous pattern of interactions on the 
cluster state. Therefore one has 

Z G (J)=Z D (J,J'). (5) 

In this way it has been shown that the 2D Ising model with complex inhomogeneous couplings is 
complete 0. 

Like any completeness result, an interesting question is whether there are other complete models. 
The situation is reminiscent of results on NP completeness of certain problems in computer science 
ll24l . In this direction it has been shown in j6l|TJ that the four dimensional Z2 lattice gauge theory, 
with real couplings, is complete for producing any spin model in any dimension. This result was 
extended in [H to show that the four dimensional U(l) lattice gauge theory with real couplings can 
produce, to arbitrary precision, a large number of continuous models (those whose Hamiltonian allow 
a finite Fourier series). Certainly there may be many other complete models which can be converted 
to each other. Exploration of the set of complete models certainly will add to our insight and to our 
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power in connecting different models with each other. 

In this paper we show that the discrete form of 4 field theory, on a two dimensional rectangu- 
lar lattice is also complete in the sense that the partition function of any continuous model on any 
graph with any type of interaction can be obtained, to arbitrary precision, from a </> 4 model with in- 
homogeneous complex couplings on an enlarged 2D lattice. 

The structure of this paper is as follows: First we gather the necessary ingredients for our analysis 
in section ©, i.e. elementary facts about Continuous Variable (CV) states, operators and measure- 
ments. We then reconsider in a new language three types of continuous variable stabilizer states, 
namely the CV Kitaev states, the CV extended Kitaev states and the ordinary graph states and the 
relations between them. In section (01 we introduce the quantum formalism for scalar field models on 
arbitrary graphs and investigate properties of these models, properties which are made transparent by 
using the quantum formalism and are otherwise not easy to unravel. Then in section (O, we show that 
the free field theory on two dimensional rectangular lattice is complete for free theories in the sense 
that from its partition function, every other free field theory on any graph can be obtained as a special 
case. Finally we show that the </> 4 field theory on 2D rectangular lattice is complete, in the sense that 
its partition function reduces to the partition function of any interacting model on any graph. As in 
the Ising case, the price that one pays is that the coupling constants of the complete model should be 
inhomgenuous and complex. 



2 Preliminaries: 

In this section we collect the preliminary materials necessary for generalization of the quantum for- 
malism to the Continuous Variable (C V) case 11251 l26l l27l l28l . First we review the definition of 
Heisenberg-Weyl group, and the way a unitary operator can be performed on CV state (a qumode) by 
measurements of an appropriate graph state. We end this section with a note on decomposition of CV 
unitary operators. 

2.1 The Heisenberg-Weyl Group 

The definition of CV stabilizer states starts with generalization of the Pauli group to the continuous 
setting [29, 3 0] [31] [32] [33). For one qumode, (a term which replaces qubit in the continuous setting) 
the resulting group is called Heisenberg-Weyl group W, whose algebra of generators is spanned by 
the coordinate and momentum operators satisfying [Q,P] = il- Thus modulo U(l) phases, the group 
W is the group of unitary operators of the form w(t, s) = e lt Q+' lsP , Since the two unitary operators 

X(s)=e- lsp and Z(t)=e lt ®, (6) 

have the simple relation 

Z(t)X(s) = e lst X(s)Z(t), (7) 

multiplication of any two elements of W can be recast in the form X(t)Z(s) modulo a phase. The 
Heisenbergy-Weyl group can be represented on the Hilbert space of one particle, spanned by the basis 
states \y) q (eigenstates of Q) or \y) p (eigenstates of P). On these states, the operators X and Z act as 
follows: 

Z{t)\y) q = e«y\y) q , X(s)\y) p = e^y),, 

Z(t)\y) p = \y + t) pi X{s)\y) q = \y + s) q . (8) 
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Remark: Hereafter we denote the states \y) q simply as \y) and use \y) p for eigenstates of P as above. 
The CV Hadamard operator is a unitary non-Hermitian operator defined as 

H := / \y) p (y\ q dy = -±= / <? xy \x)(y\dxdy, (9) 



/2tt 

from which we obtain 

HQH- 1 = P, H~ r PH = Q. (10) 

Moreover from (0 we find 

H 2 := { dx\ -x)(x\, (11) 



which leads to H 4 = I. Therefore the square Hadamard operator acts as the parity operator. This 
means that 

H 2 QH-' 2 = -Q, R-' 2 PH 2 = —P. (12) 

The n— mode Heisenberg-Weyl group is the tensor product of n copies of W, i.e. W n := W® n and 
all the above properties are naturally and straightforwardly extended to n modes. Of particular interest 
is the continuous variable CZ operator which is defined as CZ\x, y) = e lxy \x, y), with the operator 
expression 

CZ := e 1 ®®®, (13) 

and satisfying the relation 

(X(t) ® I)CZ = CZ(X(t) ® Z(t)). (14) 

We will also need a more general operator CZ, namely CZ(s) = e ls ®®® which has the following 
relation 

(X(t) <g) I)CZ{s) = CZ(s)(X(t) ® Z{st)). (15) 



2.1.1 Single mode unitary operators induced by measurements 

In this subsection, we review how a CV unitary operator can be induced on a mode by measuring a 
suitable graph state. Here we restrict ourselves to operators diagonal in the coordinate basis, since 
this is the only type of operators which we encounter in our analysis. We also use the basic result of 
MQC that certain states are complete, in the sense that by a suitable sequence of adaptive single-site 
measurements on them, any other state can be reached. 

Consider a very simple two-mode graph state as shown in figure (Q~|). The first mode is in an 
arbitrary state \<f>) = J dx<f){x)\x), the second mode is in the state |0) p := J dy\y), the two 
modes have been joined by a CZ operator (shown by a line in figure©) and so the two mode state is 

|*)ia = (<7Z)ia|0)i|O p )a = J dxdye lxy 0(x)\x,y} lt 2, (16) 

where the indices 1 and 2 refer to the modes from left to right in figure (Q]). Now we project the first 
mode on the zero momentum state (0| p l33l [341 135 36 37][38][39]. The state of the second mode will 
be 

\<f>'h -i (0 p |*)ia = -L / dye^ ( f>(x)\y) 2 = H\<j>). (17) 

V Z7T J 

Thus projection of the first mode onto a zero momentum state is equivalent to the action of the 
Hadamard operator on the state \<j>) and putting it on the second mode. This is shown in figure (Q]), 
where projection is depicted by a downward arrow and the result is depicted by an upward arrow. 
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Hg(Q)|(|>) 



Figure 1 : (Color online) Projection of the left mode on the zero momentum state is equivalent to the 
action of Hadamard operator on the right mode. Projection on the state (g\ is equivalent to the action 
of the operator Hg(Q) on the right mode. A downward arrow means projection on a state, an upward 
arrow means the resulting state. 







g(Q)l<t>> 



I*) 

Figure 2: (Color online)The measurement pattern which enacts the operator g(Q) on the right-most 
mode. 

Suppose now that we project mode 1 onto the state (g\ := — h= J dyg(y)(y\. If we note that 



(g\ = -L / dyg(y)(y\ = {0 P \g(Q), (18) 
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and note that g(Qi) commutes with [CZ)\i, we find that projecting the first mode on the state (g\ 
is equivalent to the action of the operator Hg(Q) on the state \<j>) and putting it on the second mode. 
This is shown in figure ([TJ. We can write this symbolically as 

P —>H, P g ^ Hg(Q), (19) 

where in the left hand side we show the projections and in the right hand side we show the resulting 
action on the state. In order to enact the operator g(Q), i.e. remove H from Hg(Q), we need to enact 
the operator H three times. Thus using the symbols in equation ( fT9] i, we have 

PoPoPoPg — ► H 3 (Hg(Q)) = g(Q), (20) 

which is shown in figure (ffji. 



We are now in a position to state a basic theorem [40 41 . 42 43 1 in this section. 

Thereom: Let V(Q) be a polynomial of Q with real coefficients and \<$>) be an arbitrary state of an 
appropriate chain of a cluster state. Then by projecting the modes of this cluster state on the following 
three types of states, 

(Mt)\ := / dye-**«(y\, ((3 2 (t)\ := / dye^ (y\, (&(t)| := / dye^ {y\, (21) 



we can enact any operator of the form e on the state |</>), to any desired precision. The state 

will appear on the un-projected modes of the chain. 
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Proof: First we note that enacting the operator H comes for free by projecting on the zero mo- 
mentum state (0| p . Second we use the operator identity 

e tA e tB e -tA e -tB ^ e e[A,B]+o(S) (22) 

From projection on the states we find that the operators e~ lt ®, e~ lt Q~ and e~ lt( ^ 4 can be 

obtained. In view of the existence of the Hadamard operator, the algebra of anti-Hermitian operators 
in the exponential is generated by the set {iQ, iQ 2 , iQ 4 , iP, iP 2 , iP 4 }. It is now easy to see that this 
algebra contains all monomials of Q. To show this we first note that 

[iP,iQ 4 ] ~iQ 3 , (23) 

where ~ means that we have ignored numerical factors. We then note that 

[iP, [iP 2 , iQ 4 ]] ~ i(PQ 2 + Q 2 P). (24) 

The latter operator now acts as a raising operator for powers of monomials, since 

[i(PQ 2 + Q 2 P), iQ n ]] ~ iQ n+ \ (25) 

which completes the proof. In this way we can generate any unitary operator of the form e~ %v ^', 
where V is a real polynomial of Q. 



3 Three Classes of Continuous Variable States 

Let G = (V, E) be a graph, where V and E respectively denote the set of vertices and edges. The 
graph is supposed to admit an orientation. In other words, G is the triangulation of an orientable 
manifold. This means that all the simplexes of G inherit the orientation of the original manifold in a 
consistent way. The number of vertices and edges are respectively given by | V\ and \E\. 

In this section we define three closely related continuous variable states pertaining to a given 
graph G, which we call the Kitaev state \Kc) l44l . the extended Kitaev state \Kq) and the graph 
state \^g) B51 l46l . We will then determine the mutual relationships of these states, which will play 
an important role in our proof of completeness. These are the generalizations of known states in the 
qubit case, where they have been possibly named differently in other works. For example in the qubit 
case, extended Kitaev states have been called pseudo graph states [2j, however in view of their ex- 
plicit construction and stabilizers, we think that the name Kitaev or extended Kitaev states are more 
appropriate for them. 

The crucial difference between the continuous and qubit case is the fact that the operators X(t), 
Z(t), CZ and H are not equal to their inverses. Therefore a consistent and unambiguous description 
of these states on a graph, requires that the graphs be decorated with weights and/or orientations. 
We emphasize the difference between orientation, which is a Z-2 variable and weight which is a real 
variable. We will meet the necessity of each as we go along in our definitions. 

3.1 Kitaev States 

Consider an oriented graph G = (V,E,a), where a means that arbitrary orientations have been 
assigned to the edges. Any collection of arbitrary orientations on the edges is called a decoration of 
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Figure 3: (Color online) The graph for which the Kitaev state is given in d29l >. 



the graph. We assume that modes live only on the edges E and there are no modes on the vertices V 
of this graph. The CV Kitaev state is then defined 



\K G )= / d4> x dfa ...d0jv \<k - 4>j ) , (26) 
J 

where is the edge which goes from the vertex i to the vertex j. 

It is easily verified that this state is stabilized by the following set of operators: for each vertex 
i <G V, we have 

A v (t) := [] Xf{t), (27) 

e6_E„ 

where E v denotes the set of edges incident on the vertex v and the — and + signs are used for edges 
going into and out of a vertex respectively. The reason that A v stabilizes the state (l26i l is that it simply 
shifts the variable 0„ which will be neutralized under the integration. Also for each face of the graph, 
we have 

Bf(s) := IJ Z e(s), (28) 

eedf 

where df denotes the set of edges in the boundary of / and the + and — signs are used respectively 
when the orientation of a link is equal or opposite to us when we traverse a face in the counter- 
clockwise sense. Note that traversing all the plaquettes in this sense is meaningful for an orientable 
traingulation. Here also the effect of B p (s) on the state inside the integral is to multiply it by a unit 
factor since the phases acquired by all the edges add up to zero for a closed loop. 

As an example, we have for the graph Q in figure ©, 

\ K q) = / Dx\(f>i — 02,02 - 03,04 — 03,04 - 01,04 — <f>2)a,b,c,d,e, (29) 



where the subscripts a to e, determine the position of modes in the state. The stabilizers of this state 
are then given by (127128b as follows: 

A\ := X d 1 X a , A 2 := X a 1 X~ 1 Xb, A 3 := X b 1 X~ 1 , A± := X c X^X e (30) 



and 



Bi := Z e Z a x Z d , Bn := Z c Z b 1 Z e l . (31) 



Finally we note that all Kitaev states on a given graph, corresponding to different decorations are 
related to each other by local unitary actions. In fact if we switch the arbitrary orientation on a link 
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etj, it means that the term | • • • 0j — <f>j, ■ ■ ■) in (l26l l changes to | • • • (j)j — <f>i, ■ ■ •) where the remaining 
parts of the state will remain intact. In view of (fTTT i. this switching is achieved by a local action of the 
square Hadamard operator H 2 on the edge e. We can thus write 

\Ka*')=\ <S> H e)\ K G")- (32) 

\e:a(e)^a' (e) J 

Therefore all the Kitaev states with different decorations belong to the same class of states modulo 
local actions of square Hadamard operations. 

In view of the shift invariance (f>i — > <f>i + i], the Kitaev state has a hidden multiplicative factor 
which is in fact infinite. This symmetry can be removed by fixing a gauge (in the discrete case this is a 
finite factor which causes no problem). Therefore we will define the gauge-fixed Kitaev state, denoted 
by \K G ) as follows: 

\K%) = y# 1 # 2 ...#jv^i + 02 + --^jv)0|^-^). (33) 

This gauge-fixed Kitaev state has still the same set of stabilizers. The same reasoning for the Kitaev 
state \Kq) also works here. 



3.2 Extended Kitaev States 

Let G = (V, E, a) be defined as in previous subsection. Now in addition to the modes on the edges, 
there are also modes on vertices, as in figure ©. On such a graph, there are thus two different sets 
of vertices, which we denote by V (the ones on the nodes) and Ve (the ones on the edges). Then the 
state \Kg) is given by 

\K G ) = f d<t> x d<t>2...d<j> N (g) \<t>i - 4> 3 ) (g) \(f>i}, (34) 

where is the edge which goes from vertex i to j. As an example for the graph in figure (0), we 
have 

\Kg) = J D<j>\(pi - 4>2,4>2 - 4>3,4>A - 03, 04 - 01, 04 - 4>2}a,b,c,d,e ® 101, 02, 03, 04)l,2,3,4- (35) 

The stabilizers of this extended Kitaev state are completely different from the simple Kitaev state. In 
fact they are: For each vertex v E V, 

C v :=X V 11 Xf, (36) 

where the convention for the ± signs is the same as in Kitaev state and for each edge e G £ which 
goes from v\ to V2, 

De-Z^Z^Z- 1 , (37) 

where again we have suppressed the continuous arguments of these stabilizers for ease of notation. 
For the example given in figure (|4j, some of the stabilizers are: 

C\ = X\X a X d 1 , C<i = X2X a 1 X e 1 Xi J , D a = Z\Z a 1 Z 2 , Db = Z2Z b 1 Z%, ■ ■ ■ (38) 
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Figure 4: (Color online) The graph Go for which the extended Kitaev state is as in ( 1351 1. 

3.3 Weighted Graph States 

Finally we come to the definition of continuous variable weighted graph states. Here as in the qubit 
case we start with an initial product state of the form \ = \0)f V - However there is an important 
difference in that on each edge, instead of the simple CZ operator, we can act by the CZ(J) operator 
where the real parameter J may depend on the edge. Therefore we obtain what we call a weighted 
graph state. Denoting the collection of all weights by J, we have 

I* g(J)) = <g)(CZ(J e ))\n), (39) 
where \ = \Q)f V - The explicit form of the state will then be given by 

|*gG7)} = J D<t*£«.i> J ^|0 1; ---<M, (40) 

where N is the number of vertices and < i, j > denotes the edge connecting the vertices i and j 
carrying weight J^ . Note that in contrast to a decorated edge which is denoted by (going from i 
to j), a weighted edge is denoted by the symmetric symbol (i, j) . According to (TOT l. the stabilizers of 
this state will be of the form 

Ki(t) := Xi(t) Yl Z^t), Viey, (41) 

where in the left hand side we have suppressed the dependence on the weights for simplicity. 

We emphasize that definition of the Kitaev states and extended Kitaev states require the underlying 
graph to be decorated, while a graph state needs only a weighted graph for its unambiguous definition. 
We are now left with an important question of whether there is a simple relation between the above 
three kinds of states or not. The answer turns out to be positive and is explained in the next subsection. 

3.4 Relations between the above three states 

Consider the extended Kitaev state corresponding to a decorated graph G = (V, E, a). The explicit 
form of the state is shown in ( f34b . From ([34-b it is clear that if we project all the vertices in V on the 
zero-momentum basis, we will arrive at the Kitaev state for the same graph. More explicitly we have 

<n^> = -^m\Kc~), (42) 



9 



where |fi) = lO)® 17 , and we have explicitly indicated the decoration a. 

It is instructive to understand this in an alternative way, that is by showing that measurement in the 
momentum basis actually transforms the stabilizer set of the extended Kitaev state, i.e. S(\K G )) to the 
stabilizer set of the Kitaev state, S(\K G )). From the stabilizer formalism we know that measurement 
of a state 1^) in the basis of an operator M, removes all the operators which do not commute with M 
from the set S'd*)) and leaves us with a smaller subset. This subset is generated by all the original 
generators or their products thereof which commute with M. Having this in mind, it is straightfor- 
ward to see that measurement in the momentum (the X basis) leaves all the vertex stabilizers C v intact 
(except of course removing the vertex X v from it), hence changing it to A v as in d2Tb . However since 
X v does not commute with Z v , measurements of all the vertices remove all the generators D e . The 
only combinations which survive this elimination will be their product around any faces. These are 
nothing but the operators Bf for all faces, which are just the right stabilizers of \Kg°)- 

Let us now study the relation between the extended Kitaev states and graph states. It turns out that 
there is a simple relation between the two only if the weights of the edges incident on each vertex add 
up to zero, that is, if 

J ij = °> v *- (43) 

3 

Since Jy = Jji, this means also that = 0. In such a case we can convert an extended Kitaev 

state to a weighted graph state on the same graph by suitable measurements on the edges. To this end 
we proceed as follows: Let us project each edge of the extended Kitaev state on the state (02{Uj)\ 
defined in d2TT >. 

Let (£a(t)| := U eij (A»(ttf)|, then we have 

(fo(t)\K G )= /"z^^-^-^V,"-^}- (44) 
If we choose the parameters tij of quadratures so that J^i Uj = J2j = 0' we ^ n( ^ 

(/3 2 (t)|T? G ) = J D^<^> lt ^\^,--^ N ) =: |* G (it)), (45) 

which is a weighted graph state with weights ity assigned to each edge (i, j) . 

Note that from the Kitaev state for the rectangular lattice | K\j ) , the Kitaev state for any other graph 
\Kq) can be obtained simply by measurement of the edge modes in the momentum or coordinate 
bases. In fact projecting an edge mode on the zero-momentum state |0) p removes that link from the 
graph while projecting it on the zero coordinate state \0} q merges the two endpoints of that edge. 
These are shown in figure ((5]) and are proved as follows: Let 

iv) = f widxhD'ttfa - 2 > n i^i - 4n) n \<t» - ■ ■ ■ > 

where • ■ ■ denotes all the edges which involve neither the vertex 1 nor 2. Projecting this state on the 
state p (Q\ on the edge (1,2), leaves us with 

|^p) = f d^dfcD'cj) 11 \^ - l[ |0 2 - <j> 3 ) ■ ■ ■ , (46) 
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Figure 5: (Color online) The effect of measurements of an edge in an extended Kitaev state. Measuring 
in the momentum (X) basis (i.e. projecting onto the p (0|) removes the edge, while measurements in 
the coordinate (Z) basis (i.e. projecting onto the p (0|), merges the two end points of that edge. 




Figure 6: (Color online) The Kitaev state on the rectangular lattice, when specific edges are measured 
in the X or Z bases, will turn into the Kitaev state on the Hexagonal or Triangular lattices. 



which is nothing but the same Kitaev state with the edge (1,2) totally removed. On the other hand 
projecting \ip) on the state g (0| on the edge (1, 2), leaves us with 

ieL j£R 

which means that the two endpoints of the edge (1,2) have been merged together. With these two 
simple rules of deleting and merging one can obtain the Kitaev state of any graph starting from the one 
on the rectangular lattice. Figure (O shows an important example in which measurements of some 
of the edges in the momentum basis (and hence removing them), transforms | A'n) to the Kitaev state 
on the hexagonal lattice \Kh). Measurement of the same edges in the coordinate basis (and hence 
merging the two endpoints) produces a uniform lattice whose faces are triangles, hence a triangular 
lattice. This is in accord with the fact that the Hexagonal and Triangular lattice are dual to each other, 
a subject which will be explored further in the sequel. 
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Finally we use the well-known universality of cluster states proved in the context of measurement- 
based quantum computation [[33] [34j [35] [36l [37l [38] [39] to state that both the Kitaev state \K G ) and 
the extended Kitaev state \Kg) can be obtained by Gaussian measurements from a sufficiently large 
cluster state |^n). The fact that the measurements need only be Gaussian is due to the fact that all 
these kinds of states are stabilized by subgroups of W n and hence they can be converted to each other 
by unitary operators belonging to the Clifford group. Using the well-known fact from MQC, that 
Clifford operators can be implemented by Gaussian measurements, we arrive at the proof of the above 
statement. 



4 Quantum formalism for the partition functions 

In this section we show how the partition function of a classical model defined by a Hamiltonian on a 
continuous variables on an arbitrary graph, can be expressed in the quantum formalism. First consider 
the case where there is no local term or on-site interaction, that is the Hamiltonian is of the form 

H = Y,V i} (<k-<l>i), (48) 

where Vij (x) is an arbitrary function. We allow for the function Vij to depend on the edge in order 
to cover also the inhomogeneous cases. The partition function of this model is, 

Z(G, {V}) = J D^e-^CJ) (49 ) 

where N is number of vertices and we have absorbed the parameter (3 = — ^ in the Hamiltonian. We 
will do this in all expressions of partition functions which follow. 

Remark: We have defined the partition function in the form d49) . in order to be able to deal with 
unitary operators in the measurement based quantum computation. Dealing with non-unitary opera- 
tors does not pose any problem in the Ising model Q, since states like \a) = e~° |0) + e° \1) are 
normalizable states. In the continuous case the analogue of the above state may be non-normalizable, 
rendering the projection to such states problematic. Instead we resort to partition functions of the type 
(fTob with the understanding that the final results, like dualities, and completeness can be analytically 
continued to the whole complex plane. 

Due to the shift invariance of the Hamiltonian fa — > fa + £, the above partition function is 
infinite, so we have to modify the partition function by fixing a gauge, which we will do later on. For 
the present we deal with the above partition function as it is. Defining the product state 

|a) = (g)K->, (50) 

Gij 

where 

\ aij ) = J dxe- iV «M\x) 
is defined on the edge e.y, one can then write the partition function d49b in the quantum formalism as 

Z'(G,{V}) = (a\K G ), (51) 

where \Kq) is the Kitaev state on the graph G (|26l >. In this way, as in the qubit case, the pattern of 
interactions is encoded in the entangled Kitaev state and the strength of interactions (including the 
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temperature) are encoded into the product state (a\. To fix the shift invariance we define a gauge-fixed 
partition function as 

Z(G, {V}) = f D<f>S{J2 &)e-* E <«.i> v **l+*-+i). (52) 

J i 

Note that other forms of gauge fixing terms are possible, but we will deal with this simple one. Also 
note that the shift invariance is also present in discrete models, however in those cases the multi- 
plicative factor is finite and not divergent, hence gauge fixing is not necessary. Using the gauge-fixed 
Kitaev state, we can write this in the quantum formalism as 

Z(G,{V}) = (a\K G ). (53) 



Let us now consider an edge e and insert the operator Q e inside the inner product (1531 . In view of 
(f50T > and (l53l we will have 



(54) 



G I 

where by ( ), we mean the statistical thermal average. Similarly by acting the momentum operator on 
the state (a\, we find, 

(a\Pe\K° G ) 

( a | A -0) = - <h)h ( 55 ) 

where by ' we mean derivative with respect to the argument. We will later see an important application 
of these equations when they are combined with the topological properties of the Kitaev states. 

Consider now the case where there are on-site interactions, then the Hamiltonian will be 

H = Vij (<f>i -4>i) + Ys W i(M> (56) 

and the above formalism will be extended as follows: 

Z(G,{V},{W}) = (a\K G ), (57) 

where 

l«> =®K->®I Q! <>. ( 58 > 

in which 

|ai) = / dxe- iWi &\x) 



5 Applications of the quantum formalism 

Let us now try to understand some of the properties of a continuous variable statistical model through 
the quantum formalism. Certainly the results that we will find, like duality, can also be derived by 
other means, without resorting to the quantum formalism, however this scheme make these properties 
and their derivation much more transparent. 
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Figure 7: (Color online) For any closed loop, which is the boundary of a region, equation holds for 
the correlation functions. 

5.1 Correlation functions 

Consider a graph G = (V, E) and a Hamiltonian Hq defined on it without local potentials, Wi (<fi) = 0. 
Let W~ be a closed loop on the dual graph, which is homologically trivial, (figure (|7]i) i.e. it is the 

boundary of an area. Since W~ can be written as a product of operators A s , for s inside the loop C, 
we have 

W%\K° G ) = \K° G ), (59) 

or equivalently 

Y J Pe\K° G )=0. (60) 

eec 

Then in view of OBI , this means that 

x;w&-&)>=°- (6D 

This is a general non-trivial relation which is valid for any kind of interaction V and without using 
the quantum formalism, it would have been difficult to obtain it. Note that the other kind of loop 
operator, W c defined as W c := Iliec where C is a loop in the graph, doesn't lead to a non-trivial 
relation since in view of d54i i, insertion of this operator into the inner product leads to the quantity 
X) e — fa) which identically vanishes. 

5.2 Duality 

Denote the dual graph by G. The vertices, edges and faces of G are in one to one correspondence 
with the faces, edges and vertices of G respectively. For an oriented graph we should also choose a 
convention for choosing the orientations. We choose the convention that for each edge e the dual e be 
such that the pair (e, e) form a right handed frame, as in figure dH). In view of the form of stabilizers 
of the Kitaev states (127128b and the relations d 1 0b . and the normalization of the state \K G ), we see that 

\K G ) = (2^H^ E \K G ). (62) 

Note that contrary to the qubit case the duality relation is not an involution, that is, as shown in figure 
dH), the dual of the dual of an oriented graph is not the original graph, but the original graph with all 
the orientations reversed. This is in accord with the fact that H 2 ^ I and indeed the action of H 2 on 
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Figure 8: (Color online) The dual of a the dual of a graph is the same graph with all the orientations 
reversed. 



all edges reverses their orientations. Consider now the partition function on G, with (3H as in d48l . 
We have 



Z(G,{V}) = (a\K G ) 



(27T) 



->a|.fft |ifg) = , E| N ^ G/i 



(5I/0 



(2tt) 



(63) 



where 



where 



ff|a) = fl" y dxe- lV{x) \x) = --== J 



2?r 
1 



dxdye ixv - iV W\y) = I dye~ iV ^\y), (64) 
/ dxe ixy - iV(x) . (65) 



This gives the following duality relation 



1 



Z(G,{V}) = ——Z(G,{V}), 
(2tt)V 



where e is the Fourier transform of e 1 , as given in ( I65t . An example of interest is when 

which leads to the following duality relation 

Z{G,{hj}) = - \ m * Z(G,{-^}). 



(66) 



(67) 



(68) 



6 Completeness of two dimensional 4 model for all discrete scalar 
field theories 

In this section we use the quantum formalism to show that the discrete form of two dimensional 
4 field theory is complete. Before proceeding we should make precise the meaning of the above 
statement. By the two dimensional discrete 4 theory, we mean the following Hamiltonian on a two 
dimensional square lattice with periodic boundary conditions 



H c = ^2 Kr iS ((f> r 

(r,s) 



s ) 2 + 2J h ?4r + m r <j> 2 r + <3v</> 4 , 



(69) 
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where K r . s 6 {i, —i}, i = \f— T and the real parameters h r , m r and q r denote respectively the 
inhomogeneous external field, the quadratic (mass term) and quartic coupling strengths. The linear 
terms {h r } are also necessary for completeness. Denote the partition function of this model by 

Z{U, {h}, {m}, {q}) := J D<l>e- iH % (70) 

where □ is a 2D square lattice of size N, and {h}, {m} and {q} denote the totality of all the in- 
homogeneous coupling strengths. By completeness we mean that the partition function of any other 
model 

H = Y^V r , s {4>r ~ 4>s) +^W r {<j> r ), (71) 

r.s r 

on any graph G = (V, E) is equal to the partition function for the Hamiltonian H c on the 2D square 
lattice, for some specifically chosen sets {h}, {to} and {q}. 

Note that the interaction terms in the Hamiltonian H are not necessarily of nearest neighbor type. 
In fact as we will show in the proof, H can even contain fc— body interactions, as long as k is bounded 
above by a finite constant independent of \V\. 

Consider a general Hamiltonian of the form (TTTT >. The partition function of this model can be 
written in the quantum formalism as 

Z(G,{V},{W}) = (a\K G ), (72) 

where 

(a| = 0(^10(^1, (73) 

e i,3 i 

and 

( aij \ = J 'dye- iV »M(y\, (oi| = J ' dye~ m ^ (y\. (74) 

Note that (a^l lives on the edge vertex v(eij) and (a.i\ lives on the vertex i. Rewriting the above 
states in the form 

(ay | = J dy(y\e- lV ^ = (Oje"^), (oh\ = J dy{y\e~ m ^ - (0 p \e~ iw ^K (75) 

where (0 p | is the zero momentum eignestate, we find an equivalent form for the partition function, 
namely 

Z(G, {V}, {W}) = ( p | (g) e~ w ^ (g) er m '^ \K G ), (76) 

eij i 

where (0 p | = (^) e ,. eE (0 p \ ® iG y (0 p | is the product of all zero-momentum states on the edge vertices 
and ordinary vertices of the graph G. We now note that according to d2"Tl ), we can approximate the 
unitary operators e~ lVij ^ and e~' lWi ^ to any degree of precision by a product of the operators 
H (the Hadamard) , e _it<3 , e~ lt ® 2 and e _itQ4 . As explained in ( I2.1.U , implementation of these 
operators on a state is effected by suitable possibly non-Gaussian measurements (i.e. projections 
of vertex modes on the states (f3i(t)\(i — 1,2,4)) of an appropriate enlargement of that state, i.e. 
one simply adds necessary nodes, glues them by CZ operators and measures the additional nodes as 
exemplified in figure (ffj) to affect a desired unitary gate on the original qumode of the lattice. Let us 
denote this intermediate graph by G', its associated state by I^g')' anc l tne collection of all necessary 
measurements on it by (/3i,2,4 1, then we will have 

0e-<^<3)0 e -W<S)|jr G ) = (/3 li2)4 |* G ,), (77) 

eij i 
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The configuration of the graph G' may be complicated, but the important point is that \^>g') is noth- 
ing but a stabilizer state and hence in principle it can be obtained from a cluster state by Gaussian 
measurements (f3i(t)\(i = 1,2). Therefore we have 

- (A,a|*p>. (78) 

Note that up to now all the projections have been performed on the vertices of the cluster state. The 
cluster state |^n) can be a weighted cluster state where the weights of all edges are ±1 and for each 
vertex the weights of all edges add up to 0. Such a cluster state can be obtained from an extended 
Kitaev state on the rectangular lattice by projecting all the edge- vertices on the states (^(±1)1 = 
J dye TW (y\ according to whether the weights of the edges are +1 or —1. 

This is the only place where projections are made on the edge vertices, and in fixed directions (i.e. 
eigenstates of XZ ±l ). We write this symbolically in the form 

|*p) = (±\K D ). (79) 

Combining d77l i. d78l and d79l we finally arrive at 

Z(G,{V},{W}) = (0 p |</? li2) 4,± e |^n>, (80) 

where ± e encapsulates all the measurements XZ^ which are performed on the edges of the extended 
Kitaev state, and (/?x,2,4| represents all the projections (/3j| for i = 1, 2, 4 on the vertices. Putting all 
this together we finally arrive at the result that 

Z(G, {V}, {W}) = Z(n, {h}, {m}, {q}). (81) 

It is important to note that since the edge vertices are measured in the ZX ± bases and the resulting 
edge is projected on the eigenstates |±)| oc J e v \y)dy, the interactions between neighboring ver- 
tices in the complete model is restricted to be of the type ±i(0j — 0j ) 2 - I n m is wav a ^ me couplings 
in the original model have been transferred to the mass and potential terms on the vertices. 
As a byproduct this argument shows that when the original model has only quadratic couplings (i.e. it 
is free field), then there is no linear or quartic coupling in the model on the rectangular lattice which 
reduces to this model by Gaussian measurements. This means that the free field theory on the 2D 
rectangular lattice is complete and can produce any other field theory on any lattice. 
It is a simple matter to show that the 4 theory can also reproduce models with fc— body interactions. 
We know from |8|, that the 4D U(l) lattice gauge theory is complete. Therefore it is enough to show 
that 4> 4 theory can reduce to AD U(l) lattice gauge theory. The Hamiltonian of the latter model is 
given by 

i7=-^JpCOs(0 1 -0 2 -03 + 04), (82) 

p 

where p denotes a plaquette, J p denotes the coupling constant on p, and 0j's are the continuous 
variables around p. The indices 1, 2, 3 and 4 denotes the edges of p when traversed in clockwise 
direction. The point is that the partition function of such a model can again be written as a scalar 
product Z = (a\G) =, with (a\ a product state over all plaquettes, 

(a\ = ® p {a\ p , (a\ p := J dye^^ (y\ (83) 

and \G) a new stabilizer state, 

|G> = / Dx\ ■ ■ ■ , (a?i -X2-X3+ Xi) p , ■■■). (84) 
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If we now note that (a\ p can be written as (a\ p = (0| p e lJp cos ® and the latter operator can indeed be 
expanded to any desired accuracy in terms of e~ lt ® and e~ lt ® , the assertion will be proved along 
the same line as indicated above. 

Efficiency: It is shown in |lj that simulating the partition function of any Ising or Potts type model 
on an arbitrary graph can be done on the Ising model on a square lattice with only a polynomial 
overhead in the number of spins. For more general models however, an exponential overhead may 
be needed. A similar statement is true also in our case, namely for simulating the partition function 
of a model with nearest neighbor interactions on a graph with N vertices, we need a </> 4 model on a 
cluster state with P(N) vertices where P(N) is a polynomial of N. This result is a combination of 
the universality result of cluster states which by only a polynomial overhead can produce any other 
quantum state and the fact that any quantum unitary e tV ™' can be decomposed to a product of poly- 
nomial number of unitaries of the form e ltq " , with n = 1, 2, 4 to any degree of precision. This result 
is also true when each site interacts with a finite number k of its neighbors, where k is independent of 
N. For more general models, an exponential overhead in the number of sites, will be necessary like 
the case of Ising model. 



7 Discussion 

The concept of completeness of certain statistical models is very fascinating. The idea that in prin- 
ciple a single complete model, like the 2D Ising model in its rich phase structure, various phases of 
all the other models, regardless of their lattice structure, type of order statistical variables and order 
parameter, and the interactions, is an idea which needs much exploration in the future. One of the 
basic questions is that what other types of models are complete. In this paper we have shown that the 
(f> A field theory is a complete model. Like the case of Ising model J2) or the U(l) lattice gauge theory 
0, our proof is an existence proof at present. The next step in such a program will be to show how 
other specific models can be obtained from these complete models and what insight about them can 
be obtained. Like the existence proof itself, this step relies also heavily on techniques from quantum 
information theory, notably the measurement-based quantum computation. In particular it will be very 
desirable to formulate an algorithmic approach for deriving any specific model (its graph structure and 
coupling strengths) from a complete model. In this way apparently unrelated models will be linked to 
each other and the insight gained from this approach will have far reaching consequences in statistical 
mechanics, exactly solvable models and critical phenomena. 
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